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I. Introduction 


The overall objective of this project is to advance our understanding of the three- 
dimensional structure, composition, and dynamics of Io’s atmosphere and its significant changes 
as Io moves in its orbit about Jupiter and undergoes complex, coupled interactions with the 
corotating plasma torus that shape its atmosphere, ionosphere, and loss of neutral and plasma 
species. The proposal provides a timely, unique, and very promising approach to this problem. 
It combines state-of-the-art hydrodynamic and kinetic theory atmospheric models with new 
ground-based, HST, and Galileo observations which, for the first time, provide us with spatially 
resolved images for the distribution of several atomic species (O, S, newly discovered H, and 
hopefully, Na) in Io’s atmosphere and also with ten new ionospheric profiles. It is proposed to 
undertake in a systematic way the natural extension and refinement of an already established 
multi-species hydrodynamics model with photo- and gas-phase chemistry for Io’s atmosphere 
(Wong and Johnson 1996) and a hybrid fluid/kinetic model (Marconi, Dagum and Smyth 1996) 
to establish a fundamental base upon which to study Io’s atmosphere more fully and the effects 
of the magnetosphere on the atmosphere because of plasma chemistry, ion-driven gas escape, 
and electromagnetic heating processes near the satellite. The three-year research plan is 
summarized in Table 1. 


Table 1. Three Year Research Plan 


Model 

2-D Hydrodynamic 

Year 1 

Add H chemistry, ion 
& electron processes, 
and evaluate their 
effects. 

Year 2 Year 3 

Perform support model calculations for the 3-D fluid 
and the 2-D hybrid models; add ionospheric chemistry. 

2-D Hybrid 

Implement simple 
chemistry in the model. 

Add streamlines; undertake calculations to study the 
non-equilibrium departures of the hydrodynamic model 
and the escape of gases from the atmosphere. 

3-D Hydrodynamic 

Incorporate heating, 
chemistry, electron & 
ion processes evaluated 
above in the 2-D 
model. 

Test model; undertake modeling studies for O, S, H, Na 
observations and also for different locations of Io about 
Jupiter; begin integration of streamlines from E & M 
solution; add ionospheric chemistry. 

Electromagnetic 

Interaction 


Develop code to solve for Complete development; 

the electric and magnetic undertake coupling of the 

fields and the streamlines. interaction model & the 

atmospheric model. 


II. Summary of Work Performed in the Second Year 

In the second project year, research efforts have focused on (1) evaluating the effects of 
hydrogen and chlorine chemistry in the two-dimensional hydrodynamic model for Io’s 
atmosphere, and (2) developing a three-dimensional hydrodynamic model for Io’s atmosphere. 
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Progress has been much slower than expected because of an extended illness of the Principal 
Investigator that impacted most of the project year and because of staffing changes. Much of the 
research work was therefore postponed to the third project year. 


2.1 Hydrogen and Chlorine Chemistry for Io’s Atmosphere 

Exploring the impact of new chemistry for H and Cl in Io’s atmosphere was initiated this 
year. The presence of atomic hydrogen in Io’s atmosphere was recently identified 
spectroscopically by Trafton (2000) in Lyman-oc (1216 A) emission corresponding to several 
hundred Rayleighs from measurements acquired with the Goddard High Resolution 
Spectrograph (GHRS) of the Hubble Space Telescope (HST). The likely presence of atomic 
hydrogen in Io’s atmosphere was suggested earlier by Frank and Paterson (1999) from the 
detection in Galileo spacecraft PLS data of H + pickup ions near Io and by the analysis of Galileo 
spacecraft plasma wave (PWS) data near Io by Chust et al. (1999). The detection in Io’s 
atmosphere of atomic chlorine in the emission lines at 1379.5 A and 1347.2 A and of atomic 
hydrogen in Lyman-a emission was also more recently reported by Retherford et al. (2000) from 
observations acquired with the Space Imaging Spectrograph (STIS) of HST. The likely presence 
of atomic chlorine in Io’s atmosphere was expected from detections in the plasma torus of Cl + 
emission at 8579 A from groundbased observations by Kuppers and Schneider (2000) and of Cl + 
emission at 1071 A and CP* emission at 1011 A from Far Ultraviolet Spectroscopic Explorer 
(FUSE) measurements by Feldman et al. (2000). 

If volcanoes are the primary source of chlorine species in the atmospheric, the most likely 
existing parent gas molecules have been studied by Fegley and Zolotov (2000) using 
thermochemical equilibrium calculations for an O-S-Na-K-H system. If the abundance of (Na + 
K)/C1 >1, the major Cl gases are chlorides of Na and K, the major Na gases are NaCl, Na, and 
(NaCl) 2 and the major K gases are KC1, Na, and (KC1) 2 . If the abundance of (Na + K)/C1 <1, the 
major Cl gases change to Cl 2 , Cl, S 2 C1, and SC1 2 . Fegley and Zolotov suggest that the abundance 
of H-bearing gases are likely insignificant for plausible hydrogen abundances in Io’s atmosphere. 

An initial search for known reactions involving the major species in Io’s atmosphere 
(S0 2 , SO, O, S, 0 2 ) and the minor species H and Cl (and likely parent molecules and chemical 
products) has produced the list in Table 2. The list of reactions does not contain many major- 
minor species reactions. In addition, some of the species abundances or their reaction rates are 
small so that certain reactions can be ignored in developing a closed chemical system. More 
work is, however, needed to identify additional important reactions and to determine the relevant 
rates of all reactions so as to develop an appropriately closed system for the H and Cl chemistry 
in Io’s atmosphere. 
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Table 2. Chemical Reactions for H and Cl Species 


Reactions 




Rate 

Cl 2 + hv 


Cl 

+ Cl 

Jl 

H + hv 

2 


H 

+ H 

J2 

CIO + hv 


O 

+ Cl 

J3 

NaCl + hv 

-> 

Na 

+ Cl 

J4 

HC1 + hv 

-> 

H 

+ Cl 

J5 

Cl + O 

2 


Cl 

+ CIO 

kl 

Cl + Na 

2 

-> 

Cl 

+ NaCl 

k2 

cio+ SO 

2 

-> 

Cl 

+ SO 

3 

k3 

cio+ so 


Cl 

+ so 

2 

k4 

cio+ o 

-> 

Cl 

+ O 

2 

k5 

CIO+ OH 

-> 

Cl 

+ HO 

7 

k6 

Cl + OH 


HC1+ O 

kl 

Cl + H 

2 


HC1+ H 

k8 

H +C1 


HC1 + Cl 

k9 

HCl + O 


Cl 

+ OH 

klO 

HC1 + OH 


Cl 

+ h 2 o 

kll 

H +0 

2 

-> 

H 

+ OH 

kl2 

H + OH 


H 

+ h 2 o 

kl3 


If, for example, either Cl 2 or NaCl were the parent molecule for chlorine and if the hydrogen 
species were not included, a particular simple closed system for chlorine chemistry is obtained by 
including only the reactions labeled by the rates Jl, J3, J4, kl, k2, k4, and k5, since the reaction rate k3 
is many orders of magnitude smaller than k4 and the abundance of S0 2 is no more than about an order 
of magnitude larger than SO. As a first step in assessing the impact of chlorine chemistry, this simple 
closed system has been added to the S0 2 and Na chemistry already included in the two-dimensional 
multi-species hydrodynamic model for Io’s atmosphere of Wong and Smyth (2000). The complete set 
of reactions for this chemistry is summarized in Table 3. At the end of the second project year, two- 
dimensional hydrodynamic model calculations are in progress to determine the two-dimensional 
[altitude and solar zenith angle (noon-to midnight)] structure of the chlorine species and their impact on 
the sodium species, which are coupled because of reactions J10 and kl 1 in Table 3. In the calculation, 
Cl 2 is assumed to be the parent molecule for chlorine. The surface source strength of Cl 2 will be varied 
relative to the surface source strength of the sodium bearing species to explore the relative 
compositional behaviors of these two chemical families in the atmosphere. These calculations will be 
presented (Smyth et al. 2001) at the meeting “Jupiter: The Planet, Satellites and Magnetosphere” 
scheduled for June 24-30, 2001 in Boulder, Colorado. 
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(CFDLIB). CFDLIB is a computer code developed at Los Alamos National Laboratory. It is a 
collection of hydrocodes using a common data structure and a common numerical method, for 
problems ranging from single-field, incompressible flow, to multi-species, multi-field, 
compressible flow. The computer code solves numerically the equations of transient, multi- 
material, compressible fluid mechanics. This code has been installed on an in-house SGI Power 
Challenge L computer with 10 R- 10000 processors. A test case run of has been successfully 
executed. 

In order to apply the code to the Io’s atmosphere, a 3-D spherical coordinate mesh has 
been created using two blocks, which connect to each other. One block forms the east 
hemisphere and one the west hemisphere. The other two ends of the blocks are shrunken to form 
the poles. A simple version employed in initial testing of the model had 9 boxes in the latitude 
direction, 20 boxes in the longitude direction and 30 layers from the Io’s surfaces to 500 km, 
evenly divided. The number of boxes is to be increased after the testing phase is completed. Io’s 
gravity, which varies with altitude, has been added. The code was tested for hydrostatic and 
simple dynamic states of Io’s atmosphere. For the hydrostatic state, an instability developed that 
prevented the code from running to completion. We are investigating the source of the 
instability. With the loss of the staff scientist early in the fourth quarter of this second project 
year who was developing this code, this investigation will be postponed until the third project 
year. 


The gas species S0 2 , SO, S, O, 0 2 , Na, NaS, and NaO have been included in the three- 
dimensional model. These species are the same as those adopted in the two-dimensional model 
calculations for Io’s atmosphere published by Wong and Smyth (2000). Subroutines have been 
included in the three-dimensional code for adding energy sources, and the code is now ready for 
adding plasma heating, Joule heating, solar heating, and non-LTE cooling in the atmosphere. In 
addition, subroutines for the sources and losses of chemical species have been included in the 
code. A utility based upon IDL to display graphically the three-dimensional model output has 
also been developed to facilitate visualization of current diagnosis and future model results. 
Once the three-dimensional code for Io’s atmosphere is operative, it will be verified by 
comparing two-dimensional calculations with the published two-dimensional model calculations 
of Wong and Smyth (2000). 
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